Angular momentum transfer to a star by gravitational waves 
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Summary and Introduction 

Interaction of a stocliastic background of gravitational radiation with celestial systems 
changes their dynamical elements in a random manner and give rise to secular changes 
in time(Berttoti 1973, Mashhoon etal 1981, Khosroshahi and Sobouti 1997). It has been 
speculated that a study of such secular changes might serve as a possible mean of detecting 
gravitational radiation. In this spirit we study the angular momentum transfer from a ran- 
dom background of radiation either to a rotating star or to an oscillating one. The angular 
momentum transferred to such objects by a continuous plane wave is proportional to time, 
t, and by an stochastic background is proportional to t^^^. 
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1. Interaction of a rotating star with gravitational waves 

Consider a polytrope rotating with an angular frequency fl, about the z-axis. The density 
p at position r, may be written as 

p{r,b) = po{r) + 6[pi(r) + p2(r)P2(cos^)] (1) 



where b = Vt^ /uj^g^ is an expansion parameter and uiosc = \J'^T^GpiS^) is of the order of the 
angular frequency of natural oscillations of the star. Suppose that a monochromatic plane 
gravitational wave with frequency ui falls on the star in a direction characterized by azimutal 
and polar angles a and /3, respectively. Assume the radius of the star is much smaller than 
the wavelength of the incident wave. In a Fermi coordinate system, the tidal force on unit 
mass exerted by the wave is the gradient of a scalar, 

= -V$ = f: XiKM^,], (2) 

where e <i; 1 is a certain characteristic amplitude of the wave and is independent of the 
frequency uj and hij{t) = Re{Aij{u)exp{—iLjt + ik.:K)} are the spatial components of fluctu- 
ations in the space-time metric by the incident wave, hij can be taken to be transverse and 
traceless. The torque exerted on the star, measured by a corotating observer, is 

f = ^ = / r X pF-d=^r = l^^beA{u;)u' /' p2ix)x^dx{h{t)i + h^it)]} (3a) 
at Jv oU p Jo 

where p, M, R are the mean density, mass and radius of the star, respectively, 

hi{t) — ^ sin2/3sin(Qi -I- a) cosa;^ — (7sin/3cos(Qt -I- a) sina;^, (3b) 

/i2 (^) = ^ sin 2/3 cos(Qi -|- a) cos ut^ a sin sin(Oi -|- a) sin ujt^ (3c) 

and (7 = 1, and —1 for right and left circular polarization, respectively. We substitute 
eqs.(3b,c) in (3a) and integrate with respect to time. In resonance the secular change in the 
angular momentum is proportional to t. Furthermore, this change is in the xy-plane and 
manifests itself as a precession of the rotation axis. In an isotropic background of radiation. 



the mean precession amplitude is zero. Its root mean square value, however, changes as t^/'^. 
Thus, 

^Ims — ~TT — di COS (3) da\ I fdtfdw, (4a) 
LI L\ J Jo 

and finally, 

OLs ~ /' Ph2{x)x'dxrS,{n)t for t » n-\ (4b) 

p Jo 



where Sg{uj) = e'^A{uj)'^uj'^/87rG is the energy spectral density of the background radiation. 
Equation(4) indicates a random walk for the rotation axis of the star. 

2. Interaction with normal modes of a star 

The torque of a gravitational wave on a rotating star, eqs.(3), is due to deviation from the 
spherical symmetry of the star. This asymmetry may come about because of the natural 
oscillation of the star. Helioseismic waves in different spherical harmonic numbers are promi- 
nent examples of such occurrence. This is the justification for the following analysis . 
Let p(r) , p{r) and U (r) denote the density, the pressure and the gravitational potential of 
a star in hydrostatic equilibrium. Let a mass element at r undergo an infinitesimal dis- 
placement i{r,t) from its equilibrium position. It causes small changes 5p{r,t),dp{r,t) and 
SU{r,t). The linearized Euler's equation of motion is 

- p| = V((5p) + 5pWn + pV{5U) = we , (5) 
where, 

5p = -V.(pO, (6a) 



V'^{SU) = -inGSp. (6c) 
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The displacement ^ belongs to a function space Ti, in which the inner product is defined as 
iViPi) = I PV*-^ d?^ = finite, {,r/ e ?Y. The operator W is self-adjoint on H and gives 
rise to the eigenvalue problem = <^^pCn, where lu^ are real. Using a gauged version 
of Helmholtz's theorem, one may decompose a general vector into an irrotational and a 
"weighted" solenoidal component. Thus 

e = + ^9, (7a) 
where 

Cp = -Vxp; with V X = 0, (7b) 



= p-V X V X (rxg); with V.{pQ = 0. (7c) 

Here r is the unit vector in r direction, and Xp and Xg are two scalars. See Sobouti(1980) 
for details of eqs.(5-7). As in eqs.(3) the torque exerted by a gravitational wave incident on 
a star is, 
rlJ r 

^ = dprx f'^'^d^x. (8) 
dt Jv 

From eqs.(6a) and (7), 5p for is zero and for is 

= P{r){x; + (1 + ^-K - ^^^^Xp}YUe, 0)e— (9) 

Only / = 2 modes will contribute to the torque and for simplicity we will consider the m = 
case. Equation(8) reduces to an expression similar to that of eqs.(3), 

^ = -\^^A{uj)uj'' f p{x){x!; + (- + -)x; - ^Xp]x^dx X {h\{t)i + h'^it)]} cosuj, 

(lOa) 

where 

h[ (t) — ~ sin 2/3 sin q; cos cut + a sin (3 cos a sin cut ( 10b) 
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^2(0 — — - sin 2/3sin acosa;^ — (Tsin/3cosQ;sina;t (10c) 

The root mean square of the amphtude of precession induced by an isotropic background 
of radiation becomes 

- /' P(^){X" + (- + -)X' - \x}x'dx\^S,{u:^)i (11) 

for t:^uj-' 

A numerical evaluation of eq.(ll) involves the following steps. 

a) One calculates the p and g modes of the model star belonging to / = 2. A p— mode will 
have a large irrotational component as indicated in eqs(l). For a g^— mode the situation will 
be the opposite. 

b) One extracts the irrotational component, Xp each mode, substitutes in eq.(ll) and 
carries out the integration numerically. Numerical evaluations of the rms precession ampli- 
tudes are in progress. 
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